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Abstract 

A systematic comparison of different isotopic temperatures with the thermody- 
namical temperature of a multifragment system is made on the basis of the Statis- 
tical Multifragmentation Model. It is demonstrated that isotopic temperatures are 
strongly affected by the secondary decays of hot primary fragments and the popu- 
lation of particle-stable excited states in final fragments. The He-Li temperatures, 
measured recently by the ALADIN group, are reproduced fairly well both as a func- 
tion of excitation energy and bound charge. Our analysis confirms the anomaly in 
the nuclear caloric curve. 

PACS numbers: 25.70-z, 25.70Pq, 24.60-k 

Nuclear multifragmentation in heavy-ion reactions is intensively studied at present 
both theoretically and experimentally. One of the main goals is to investigate properties 
of nuclear matter away from the ground state. A most interesting question here is how 
multifragmentation is related to a liquid-gas phase transition in a finite nuclear system. 
To answer this question one needs observables which bring information about the thermo- 
dynamical state of the system, in particular, its excitation energy and temperature. Then 
a phase transition should manifest itself by an anomaly in the caloric curve, i.e. temper- 
ature as a function of excitation energy. According to the statistical model prediction , 
the nuclear caloric curve behaves like in an ordinary liquid-gas phase transition: initially 
the temperature increases, at excitation energies between 3 and 10 MeV/nucleon it stays 



almost constant at about 5-6 MeV, and then grows again. The first regime corresponds 
to the compound nucleus (liquid phase), the second one, to the multi-fragment mixture 
(coexistence phase), and the third one, to an assembly of nucleons and hghtest clusters 
(gaseous phase). 

The first measurements of the nuclear caloric curve have been made only recently by 
the ALADIN group 0. They indeed revealed an anomalous behaviour of the nuclear 
caloric curve similar to that predicted by the statistical model In the experiment 
the so-called isotopic temperature, Tisot, was extracted from the double ratio of helium 
and lithium isotope yields. At present nuclear temperature measurements are in fast 
progress. Several groups have reported results on nuclear caloric curves for different reac- 
tions and with different isotope thermometers 0, ^, ^, ^. Therefore, it is very important 
now to understand how these isotopic temperatures are related to the thermodynamical 
temperatures of excited nuclear systems at the stage of their break-up. 

According to the method suggested by Albergo et al. [0], the isotopic temperature is 
expressed through the double ratio of isotope yields as 

Tisot = 7 ^ ■ (1) 

ln(a • R) 

Herei?= (Yi/Y2)/{Y3/Y^), B = (B1-B2)- (53-^4). Yi and Bi are the i-th isotope yield 
and binding energy, a is a constant determined by spin degeneracy factors and masses of 
the isotopes. The indexes i =1, 2, 3 and 4 refer to the isotopes with masses and charges 
(A,Z), (A+1,Z), (A',Z') and (A' + 1,Z'), respectively. 

It is clear that this expression corresponds to the grand canonical approximation as- 
suming thermal and chemical equilibrium. Moreover, it is assumed that all fragments are 
produced simultaneously at the same T and only in their ground states. These assump- 
tions are too crude for finite and highly excited nuclear systems under consideration. A 
more realistic approach should include at least two important modifications: first, the 
microcanonical treatment of the break-up channels i.e. taking into account exact con- 
servation laws for baryon number, charge and energy, and second, the feeding of isotope 
yields from the de-excitation of hot primary fragments after the break-up. The impor- 
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tance of secondary decays was demonstrated earlier by several authors (see e.g. refs. 
[§, ^, |10|)- Statistical models of multifragmentation (see reviews [0, provide a nat- 
ural framework for introducing these modifications. These models are very successful in 
describing many observed characteristics of nuclear multifragmentation (see examples in 



refs. Hill). 

The calculations below are made within the standard version of the Statistical Mul- 
tifragmentation Model (SMM) which was used for the first calculation of the nuclear 
caloric curve M. Here we outline only some general features of the model (see details in 



ref. |]T2[)- It is assumed that at the break-up time the system consists of primary hot 
fragments and nucleons in thermal equilibrium. Each break-up channel, /, is specified 
by the multiplicities of different species, N^z, which are constrained by the total baryon 
number Ao and charge Zq of the system. The probabilities of different break-up chan- 
nels are calculated in an approximate microcanonical way according to their statistical 
weights, Wf oc exp [Sf{E*, V, Aq, Zq)], where Sf is the entropy of a channel / at excitation 
energy E* and break-up volume V. 

Translational degrees of freedom of individual fragments are described by the Boltz- 
mann statistics while the internal excitations of fragments with A > 4 are calculated 
within the liquid-drop model with Fermi-gas level density. An ensemble of microscopic 
states corresponding to the break-up channel / is characterized by a temperature Ty which 
is determined from the energy balance equation 

^T{m - 1) + E Eaz{T)Naz + Ef{V) -Qf = E* . (2) 

Here m = J2 N^z is the total fragment multiplicity, first term comes from the translational 
motion, second term includes internal excitation energies of individual fragments and third 
term is the Coulomb interaction energy, Q/ is the Q-value of the channel /. The excitation 
energy E* is measured with respect to the ground state of the compound nucleus {Aq.Zq). 
In our semi-microcanonical treatment E* is fixed for all fragmentation channels while the 
temperature Tf fluctuates from channel to channel. 

The total break-up volume is parametrized as \^ = (1 + k)Vo, where Vq is the com- 



pound nucleus volume at normal density and the model parameter k is the same for all 
channels. The choice of k is motivated by the requirements: a) to avoid overlaps between 
the fragments and b) to provide a sufficient reduction of the Coulomb barrier, as seen in 
the kinetic energy spectra. The entropy associated with the translational motion of frag- 
ments is determined by the "free" volume, Vf, which incorporates the excluded volume 
effects. In general Vf depends on the break-up channel and therefore cannot be fixed to a 
constant, kVo, as often assumed. In the SMM we parametrize V/(m) in such a way that 
it grows almost linearly with the primary fragment multiplicity m or, equivalently, with 
the excitation energy e* = E* /Aq of the system |T^. According to this parametrization, 
V/(m) vanishes for the compound nucleus (m = 1) and increases to about 2Vq at e* R:^ 10 
MeV/nucleon. 

At given inputs Aq, Zq and e* the individual multi-fragment configurations are gener- 
ated by the Monte Carlo method. After the break-up hot primary fragments loose their 



excitation. The most important de-excitation mechanisms included in the SMM |jT2| are 
the simultaneous Fermi break-up of lighter fragments {A < 16) and the evaporation from 
heavier fragments, including the compound-like residues. In this respect SMM essentially 
differs from the QSM type models [|1^ where the compound-like channels are completely 
ignored (see discussion in ref. [Il5|). 

Now we turn to numerical simulations of the multifragmentation on the basis of SMM. 
First of all we present results for a well defined source i.e. an excited ^^^Au nucleus. The 
caloric curve is calculated by first solving eq. (2) for each particular channel and then 
averaging Tf over a large number of break-up channels. In Fig. 1 (top) different curves cor- 
respond to different choices of volume parameters. If the standard parametrization Vf{m) 
is used, the caloric curve is quite flat in the e* region between 3 and 10 MeV/nucleon. 
This is a signature of a large heat capacity in the transition region. Even a backbending 
is possible if the total volume V is not very large, say only 3Vo (k=2). In contrast, if the 
free volume would be flxed to Vf = kVq for all channels, the temperature would increase 
gradually with e*. Nevertheless, as seen in Fig. 1, some flattening in the caloric curve 
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is predicted also in this case. The reason for the different behaviour is clear: at e* < 10 
MeV/nucleon the multiplicity-dependent free volume V/(m) is smaller than kVq that leads 
to a higher temperature of the system. In the following calculations we use k=2 and the 
standard SMM parametrization of V/(m), which gives a plateau in the caloric curve. As 
we will see below such an behaviour is favoured by the data. 

The characteristics of the system change drastically when e* increases from 3 to 10 
MeV/nucleon. In the lower part of Fig. 1 we display several observables calculated after 
the completion of all secondary decays. A heavy residue, usual for the evaporation-like 
processes, practically disappears. This is signalled by the maximum fragment charge, 
Zmaxi which drops rapidly from 60 to about 6 in this region. At the same time, the 
number of Intermediate Mass Fragments (IMFs: 3 < Z < 20) first increases and then 
goes through the maximum, Nimj ~ 8, in the end of this region. The multiplicity of all 
charged particles, Nch-, grows with e* almost linearly. It is interesting to note that in the 
transition region the number of free neutrons, Nneui is almost constant and close to the 
neutron excess in the initial ^^''Au nucleus. This happens because the system breaks up 
predominantly into fragments with N ^ Z (see also [1^). By comparing upper and lower 
parts of Fig. 1 one can conclude that the neutron multiplicity is nearly proportional to 
the temperature of the system but not to the excitation energy. 

We have also calculated final isotope yields in the disintegration of ^^^Au nucleus. 
Several isotopic temperatures were obtained by applying formula (1) to different isotope 
pairs. Results are shown in Fig. 2 together with the microcanonical temperature Tmicr- 
One can see that the plateau is almost washed out and all isotopic temperatures increase 
gradually with e* . This behaviour can be explained by the de-excitation of hot primary 
fragments leading to their cooling and side-feeding of isotope yields. Since the energy 
conservation is controlled at all stages of the calculations, the SMM leads naturally to the 
cooling of emitters in endothermic processes responsible for the fragment de-excitation. 
In the case of sequential evaporation the first fragments are emitted from a source char- 
acterized by the emission temperature Tmicr- But the next generation of fragments comes 



5 



from a cooler residue leading to a lower apparent temperature 0. This cooling mech- 
anism can explain partly the difference between the isotopic temperatures and Tmicr at 
lower excitation energies (e* = 1 ^ 6 MeV/nucleon), when heavy residues {Z^ax > 20) 
survive in the break-up. At e* > 3 MeV/nucleon another de-excitation mechanism be- 
comes increasingly important, i.e. the one-step Fermi break-up where only particle-stable 
decay products are allowed. It is mainly responsible for the production of light isotopes, 
in particular He and Li, through the deep disintegration of heavier fragments {A ~ 15). 
Since the available energy (per nucleon) is considerably lower in this process than in the 
primary break-up, the apparent isotopic temperatures are also lower. Finally, at high ex- 
citation energies, e* > 10 Mev/nucleon, when predominantly light fragments are formed, 
the reduction of the available energy for secondary break-up becomes less important and 
isotopic temperatures, in average, approach Tmicr- 

From Fig. 2 one can also see that the temperature measurements can be significantly 
obscured by the irregularities in the excited states of light fragments. In our standard 
calculations the final isotope yields include the fragments in particle-stable ground and 
excited states decaying by the 7-emission. For the considered isotopes they are: 3.56 
MeV for ^Li, 0.48 MeV for ^Li, 0.43 MeV for ^Li, 3.37, 5.96, 6.18 and 6.26 MeV for i°Be. 
No such excited states are seen in ^''^He and ^Be and therefore only ground states are 
included for these nuclei. One can see that the deviations from the true temperature are 
especially large (curve a) in the case when one of the isotopes has many and the other, 
only a few or no excited states, like e.g. in the ^'^Be-^Be pair. If the excited states in 
^°Be are artificially switched off, the corresponding isotopic temperature (curve b) changes 
drastically and follows the common trend. To suppress the fluctuations associated with 
the nuclear structure effects one can use an ensemble of isotope thermometers as suggested 
in ref. [Q. Another possibility is to use isotope pairs with only a few low-lying states 
which almost compensate each other, like e.g. ^Li-^Li. In this respect it is preferable to use 
thermometers with isotopes not heavier than lithium, such as the He-Li one. But in this 
case one is facing another problem, i.e. the contamination of yields by the pre-equilibrium 
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emission prior to the break-up. This contribution is most important for hghter fragments 
and can be evaluated only on the basis of dynamical simulations. 

To apply SMM for analyzing experimental data one needs to know the characteristics 
(masses, charges, excitation energies) of thermalized emitting sources. A clear identifi- 



cation of such sources is made only in a few cases. One example is given in ref. |14 | 
where the emitting source with mass 315, charge 126 and thermal excitation energy of 
about 5 MeV/nucleon was found for central Au+Au collisions at 35 AMeV. The SMM 
calculations reproduce nicely the fragment charge distribution for this reaction yielding 
the emission temperature of about 6 MeV|T^. By inspecting Fig. 2 one can see that 
there is no contradiction between the value T(^'''Li/^'^He)~ 4.6 MeV measured for this 
reaction and the SMM prediction of T^icr ~ 6 MeV. In accordance with experiment 
is also that the Be-Li temperature is much higher than other isotopic temperatures. On 
the other hand, two other isotopic temperatures presented in Fig. 2, T(^'^H/'^'^He) and 
T('''^Li/^'^He), are predicted too high and in inverse order compared to T(^'''Li/'^'^He) 
[^]. Our analysis shows that the yields of neutron-rich isotopes, such as '^H and '^Li, are 
quite sensitive to the N/Z ratio in the decaying thermalozed source. The results of Fig. 2 
correspond to the 197Au nucleus with N/Z ^ 1.5. The correct ordering of the isotopic 
temperatures can be achieved by adjusting the N/Z ratio in the source. 

Finally we present our analysis of the ALADIN data 0, Q for peripheral Au+Au 
collisions at 0.6 and 1 AGeV (see also 0). In these experiments only fragments from the 
projectile spectators were detected. Therefore, here we are dealing with a wide ensemble of 
emitting sources produced at different impact parameters. As known JT^, the masses and 
excitation energies of these sources are strongly affected by the pre-equilibrium emission. 
Nevertheless, the ensemble of thermalized sources can be reconstructed by backtracing 



the measured characteristics of produced fragments [|13| p!q] . 

In the SMM calculations presented in Fig. 3 we have considered two different ensembles 
of emitting sources obtained in refs. and ||13| . The "experimental" ensemble of ref. 
has a wider distribution in excitation energy (up to about e* ~ 14 MeV/nucleon) than 
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the "theoretical" ensemble of ref. |]T3[ which is limited at e* ~ 8 MeV/nucleon. As seen 



from Fig. 3 the observed He-Li temperatures are better reproduced by the experimental 
ensemble. But this ensemble is certainly contaminated by the early emitted H and He 
fragments which were not separated in the data analysis. Obviously their admixture 
is larger at higher excitation energies. On the contrary, in ref. ||13| the sources were 
reconstructed by using the characteristics of fragments with Z >3 which are less affected 
by the pre-equilibrium emission. Therefore, we expect that after separating early emitted 
H and He fragments experimental points will shift closer to the prediction of the theoretical 
ensemble. 

For both ensembles we get a more steep increase of the isotopic temperature with 
excitation energy than the experimental data show (Fig. 3, top). Also within the present 
version of SMM we cannot reproduce the low temperatures extracted from the relative 
level population in light fragments such as ^Li. Our preliminary calculations show that 
the agreement with experiment can be improved by reducing excitation energies of pri- 
mary fragments and thus suppressing their secondary decay contribution. This and other 
modifications of the model are under investigation now. 

In conclusion, on the basis of SMM we have demonstrated that the secondary de- 
excitation processes and irregularities in particle-stable excited states of fragments may 
cause significant deviations of isotopic temperatures from the thermodynamical temper- 
ature of the decaying system. Our analysis shows that the ALADIN data are consistent 
with the anomaly in the nuclear caloric curve. For future studies of the nuclear caloric 
curve it is very important to separate the contribution of light clusters emitted at early 
non-equilibrium stages of the reaction. Therefore, the determination of the tempera- 
ture and excitation energy should be accompanied by a thorough kinematical analysis of 
emitting sources and fragment spectra. 
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Figure captions 

Fig.l: Top: Caloric curves as predicted by the SMM simulations for an excited ^^'''Au 
nucleus. Results are shown for four different choices of volume parameters characterizing 
the break-up configuration (see the text). Bottom: Some observable characteristics as 

9 



functions of excitation energy in multifragmentation of ^^^Au nucleus after de-excitation 
of primary fragments. "Lmax is the maximum fragment charge, A^jm/ is the multiphcity 
of intermediate mass fragments (3 < Z < 20), Nch and Nneu are the total numbers of 
charged particles and free neutrons. 

Fig. 2: Isotopic temperatures for four isotope pairs (indicated in the figure) versus 
excitation energy calculated for ^^^Au by applying formula (1) to final isotope yields. The 
microcanonical temperature of the decaying nucleus is the solid line. 

Fig. 3: He-Li temperatures (scaled by factor 1.2) versus excitation energy e* (top) and 
bound charge Zhound (bottom) for projectile spectators produced in Au + Au collisions 
at 0.6 and 1.0 AGeV. Symbols represent the ALADIN data for 0.6 AGeV 0, |l6l (dots) 
and 1.0 AGeV (triangles). The SMM calculations are made for two ensembles of 



thermalized sources: 1 from ref. 131 and 2 from ref. 
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